Novel phase transitions in XY Antiferromagnets on Plane Triangulations 
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Using Monte Carlo simulations and finite-size scaling, we investigate the XY antiferromagnet on 
the triangular, Union Jack and bisected-hexagonal lattices, and in each case find both Ising and 
Kosterlitz-Thouless transitions. As is well-known, on the triangular lattice, as the temperature 
decreases the system develops chiral order for temperatures T < T c , and then quasi-long-range 
magnetic order on its sublattices when T < Ts, with T s < T c . The behavior T s < T c is predicted 
by theoretical arguments due to Korshunov based on the unbinding of kink-antikink pairs. 
On the Union Jack and bisected-hexagonal lattices, by contrast, we find that as T decreases the 
magnetizations on some of the sublattices become quasi-long-range ordered at a temperature T s > 
T c , before chiral order develops. In some cases, the sublattice spins then undergo a second transition, 
of Ising type, separating two quasi-long-range ordered phases. On the Union Jack lattice, the 
magnetization on the degree-4 sublattice remains disordered until Tc and then undergoes an Ising 
transition to a quasi-long-range ordered phase. 
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Two-dimensional fully-frustrated XY (FFXY) models 
have been the subject of considerable interest over the 
past three decades; for a recent review see 0] ■ In addition 
to their intrinsic theoretical importance within the field 
of critical phenomena, such models can also be realized 
experimentally via Josephson Junction Arrays (JJAs) in 
a uniform magnetic field 0, H| . 

On the triangular lattice, or indeed any triangulation 
of the plane, the FFXY model coincides with the usual 
antiferromagnetic XY (AFXY) model 0, ||| , defined by 
the reduced Hamiltonian 





H = J S i ■ S i = J S 



COS 



(1) 



where Sj — (cos 6i, sin 6i) is a planar spin with unit length 
on lattice site i, J = 1/T is the inverse temperature, and 
the summation is over all nearest-neighbor pairs of sites. 

The ground states of the triangular-lattice AFXY 
model are such that the spins on each of the three equiv- 
alent triangular sublattices are perfectly aligned, with 
adjacent spins differing by an angle ±2ir/3. In addi- 
tion to an SO (2) rotational degeneracy, generic of XY 
models, there is a discrete Z2 — 0(2)/SO(2) reflection 
degeneracy, induced by frustration, corresponding to the 
two possible chiralities of each elementary triangular face. 
The chirality of a face refers to the sign of the rotation of 
the spins as one traverses the face counterclockwise; see 
Fig. Q] At positive temperatures, the Mermin- Wagner 
theorem forbids the sublattice spins from ordering. 
However, for low temperatures, the sublattice magnetiza- 
tions exhibit quasi-long-range (QLR) order (algebraically 
decaying correlations), while the chiral degrees of free- 
dom exhibit a genuine long-range order. 

Despite some early controversy [H, 0, l^4l2T[. there is 



FIG. 1: (Color online) (Left) A ground-state configuration of 
the AFXY model on the UJ lattice. Symbols ' + ' and '-' de- 
note the chirality of the elementary triangular faces. (Right) 
A BH lattice. The tripartition of the vertex set is shown in 
purple /cyan /gray. 



now a consensus [H 0, Iui - fl9l ] that for a number of FFXY 
models sharing the same ground-state degeneracies, in- 
cluding the triangular-lattice AFXY model and square- 
lattice FFXY model (20rl22l |. the phase transition asso- 
ciated with the magnetic order parameter occurs at a 
temperature T s strictly below the transition point T c for 
the chiral order parameter. It is generally accepted Q 
that the sublattice magnetizations disorder via a stan- 
dard Kosterlitz-Thouless (KT) transition [23 - 25 [ , while 
the chiralities undergo an Ising transition. The behavior 
T s < T c is supported by theoretical arguments based on 
the unbinding of kink-antikink pairs fll l26|. 

In this Letter, we study the AFXY model on the Union 
Jack (UJ) and bisected-hexagonal (BH) lattices, using 
Monte Carlo simulations and finite-size scaling. The UJ 
and BH lattices are plane triangulations that share many 
of the properties of the triangular lattice; see Fig. [T] In 
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particular, the ground-state degeneracies of the AFXY 
model on each of these three lattices arc identical. How- 
ever, as we show, the critical behavior of the AFXY 
model on the UJ and BH lattices is qualitatively different 
to the triangular-lattice case. In particular, we find that 
on both the UJ and BH lattices T s > T c . 

Like the triangular lattice, the UJ and BH lattices are 
tripartite; each consists of three independent sublattices, 
which we label A, B, C. Since each sublattice is regular, 
we use a subscript to indicate its coordination number. 
The triangular lattice consists of sublattices A$, B§ and 
Cq, UJ consists of A4, B$, Cg, and BH consists of A4, 
Be, C12; see Fig.[T] Unlike the triangular lattice, the sub- 
lattices of the UJ and BH lattices are not all equivalent, 
and this leads to interesting new physics. 

On the UJ lattice, we find that the spins on the A4 
sublattice become QLR-ordered at a different tempera- 
ture to the spins on the B 8 and Cg sublattices. Analogous 
behavior has very recently been observed in the 4-state 
Potts antiferromagnet 27|, l28|. Perhaps more surpris- 
ingly, the magnetic transition of the A4 spins, separating 
the disordered and QLR-ordered phases, appears to be 
an Ising transition. Such novel Ising transitions lead- 
ing to QLR order, rather than genuine long-range order, 
have very recently been observed [29| in generalized XY 
models whose Hamiltonians contain ncmatic-likc interac- 
tions \MM 



To our knowledge, our study presents the 
first observation of an Ising magnetic transition in the 
standard XY antiferromagnet. 

Summary of Results. — A summary of the qualita- 
tive behaviour of each lattice is presented in Fig. [U On 
each lattice, we observe a chiral transition at a temper- 
ature T = T c , with T c (Tri) > T C (UJ) > T C (BH). In 
each case, we find strong evidence that the chiral tran- 
sition is in the Ising universality class. In addition, in 
each case we also observe magnetic spin transitions at a 
temperature T s ^ T c , with T s (Tri) < T S (UJ) < T S (BH). 
The qualitative features of the magnetic transition are 
highly lattice-dependent however. On the triangular lat- 
tice, we observe a magnetic transition from disorder to 
QLR order on each of the three equivalent sublattices 
at T s > T c . The transition is consistent with the KT 
universality class. 

On the UJ lattice, we observe that the spins on sub- 
lattices B 8 and C$ undergo a transition from disorder to 
QLR order, which is again consistent with KT behaviour. 
However on the UJ lattice, we find very clear evidence 
that T s > T c . Furthermore, we find no evidence of a 
spin transition on the A4 sublattice at T s , but we find 
strong evidence that the spins on the A4 sublattice un- 
dergo an Ising transition from disorder to QLR order at 
(or extremely close to) T c . In addition, it appears that 
the spins on the B& and C$ sublattices undergo a sec- 
ond transition at T c , which separates two QLR-ordered 
phases which have different exponents but are otherwise 
equivalent. This second transition appears to be in the 
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FIG. 2: (Color online) Summary of phases and transitions 
for AFXY model on triangular (Tri), Union Jack (UJ) and 
bisected-hexagonal (BH) lattices. Red refers to disorder, blue 
to order, and green to QLR order. Temperature increases 
from left to right. 



Ising universality class, but the measured magnetic ex- 
ponent takes a non-Ising value of 1.94(1). 

On the BH lattice, the spins on each sublattice A4, 
Bq and C12, appear to transition from disorder to QLR 
order at a common value of T s > T c . The transition 
appears to be consistent with the KT universality class. 
The spins on A4 and Bq then undergo a second transition, 
consistent with the Ising universality class, at T = T c , 
while the spins on C\i do not. The C12 sublattice appears 
to display generic KT behaviour. 

Monte Carlo Simulations. — Our simulations used a 
local algorithm based on a mixture of standard local 
Metropolis [33| updates together with overrelaxed 34, 35 1 
(microcanonical) updates. Periodic boundary conditions 
are applied, and the maximum linear system size is 
L max = 384, 256 and 256 for the triangular, UJ, and BH 
lattices, respectively. The total CPU time for the simu- 
lations was approximately 60 years, with a 3.2 GHz Xeon 
EM64T processor. 

For each sublattice S = A, B, C, we define the mag- 



netic order parameter Mj = (1/Us) J2ies Si > w ^ tn 
\S\. The chiral order parameter is defined as 



M r , = 



3V3Y, 



A 



[mn(0< -e 3 )+ wa{0 3 -°k)+ sin(0 fc - 0<)] , 



where the summation is over the Ya triangular faces. 
The sequence k) is alternately chosen clockwise and 
counterclockwise on neighboring triangles. The mag- 
netic and chiral susceptibilities are then defined to be 
Xs = V s {Mj) and \c = Y A (Af c 2 ), and we also de- 
fine dimensionless ratios Qs = (M|) 2 /(M|) and Q c = 
(Mr) 2 /{Mr). In addition, on each sublattice S, we 
measure the nearest-neighbor spin-spin correlation func- 
tion Es, and define a specific- heat-like quantity C5 = 
V S ((EI) - (E s y). 

Triangular lattice. — As T decreases, the Q c data for 
different system sizes display an approximately common 
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intersection near T ps 0.512. For T is 0.512, the Q c value 
quickly approaches 1, and (Mj) converges to a nonzero 
value as L increases, implying the occurrence of chiral 
order. The spins on the three sublattices remain disor- 
dered until T ps 0.50, where the Qs data indicate a mag- 
netic transition on each of the three equivalent sublattices 
S = A 6 ,B 6 , Cg. In the low-temperature region T J§ 0.50, 
as L increases, the Qs data converge to a line of non- 
trivial T-dependent values, implying that each sublattice 
is QLR ordered. 

For small A c = T — T c , the ratio Q c is expected to 
behave like 



= Q* + ai A c L Vt + a 2 A 2 c L 2yt + b/L, 



(2) 



where a±, a 2l and b are free parameters, and yt = \jv 
is the leading thermal renormalization exponent. We 
performed least-squares fits of (J5J) to the Q c data near 
T ps 0.512. The data in the range 24 < L < 192 are 
well-described by ©, and the fit yields T c = 0.5123(2), 
y t = 0.99(1), and Q* = 0.8587(3). 

We also performed least-squares fits of the Xc data to 



L 2yh 



(ao 



^A^ +a 2 A 2 c L 2yt +b/L) , (3) 



with the fixed Ising value yt = 1. This yields a criti- 
cal point T c = 0.5123(1) and magnetic exponent yh = 
1.874(3). 

The estimated critical exponents, yt = 0.99(1) and 
ijh = 1.874(3), agree well with the exact results 2/4 = 1 
and yh = 15/8 for the two-dimensional Ising model. The 
critical value Q* — 0.8587(3) is also consistent with the 
existing result Q* sing = 0.85872528(3)for the ferromag- 
netic triangular-lattice Ising model |36j. We conclude 
that the chiral transition is in the Ising universality class. 
We note that our estimate of y t is inconsistent with sev- 
eral earlier values reported in the literature [3, IH, 37 1 
which suggested yt ~ 1.2. Our estimate of T c is close to 
the recent estimate of T c = 0.51254(3). 

Compared with the chiral transition, it is much more 
difficult to locate the spin transition of the sublattice 
magnetizations. Fitting the ansatz ([2]) to the Qa^ 
data we obtain T s = 0.50(1). It is known [38|, |39| 
that at a KT transition the susceptibility diverges as 
X oc i 7 / 4 (lni) 1 / 8 . Assuming that the spin transi- 
tion is of KT type, we therefore consider the curves 
of x J 4 6 L _7 / 4 (lni) -1 / 8 versus T for a number of fixed 
L. From the scaling of the intersections for various 
32 < L < 256 and 0.500 < T < 0.505, we obtain 
T s = 0.5040(3). This value is consistent with the recent 
estimate [H of T s = 0.504(1). 

Union- Jack lattice. — From Fig. [3l which plots 
XB a L^ 7 ^ 4 and Qb 8 vs T, we see that the spins on the 
sublattice B$ (and also C$) undergo two distinct phase 
transitions; one at T ss 0.43 and another at T ps 0.64. 
By contrast, Fig. 2] shows that the spins on sublattice A4 
undergo a single transition near T ss 0.43. 




FIG. 3: (Color online) Quantity Xb 8 L 7,/4 vs T for various 
L, on the UJ lattice. The inset shows Qb s vs T. 




FIG. 4: (Color online) Quantity xa^L 7 ^ 4 vs T for various 
L, on the UJ lattice. The inset shows XB g L~ 91 ^ 4S vs T. 



Assuming that the spin transition near T fa 0.64 is a 
KT transition, the xb 8 data were analyzed in an analo- 
gous manner to that described for the xa b data on the tri- 
angular lattice. This yields an estimate of the spin tran- 
sition point on the B$ and Cs lattices of T s = 0.639(2). 

For the chiral transition, fitting the Q c data to ([2]) 
analogously to the triangular case yields T c = 0.4316(1), 
y t = 0.99(1), and Q* = 0.8562(3). Similarly, fitting X c 
to © with y t = l fixed yields T c = 0.4316(1) and y h = 
1.874(3). This again suggests that the chiral transition 
is in the Ising universality class. 

Perhaps surprisingly, fitting the Qa 4 and xa 4 data 
near T ps 0.64 in a similar way yields Ta 4 = 0.4316(1), 
y t = 1.00(1) and y h = 1.874(2). This suggests two 
things. Firstly, within the resolution of our simulations, 
the transition points of the A4 and chiral transitions co- 
incide, implying that either T^ 4 = T c exactly, or that 
\Ta 4 — T c \ ;$ 0.0001. Secondly, the critical behaviour as- 
sociated with the A4 spin transition is consistent with the 
Ising universality class. This suggests that the magneti- 
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FIG. 5: (Color online) Cb % versus T on the UJ lattice. The 
inset shows C'a 4 - 
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FIG. 6: (Color online) Log-log plot of Xb$L~ 2 versus L^ 1 on 
the UJ lattice. The dotted lines correspond to the critical 
points. 



zation on the A4 sublattice undergoes an Ising transition 
from a disordered phase to a QLR-ordercd phase at (or 
very close to) T = T c . Based on these observations, we 
conjecture that the A4 sublattice spins undergo an Ising 
transition precisely at the chiral transition temperature 
T 

1 c ■ 

In addition, we find that the spins on sublattices Bg 
and Cg also undergo a transition near (or at) T c , in addi- 
tion to the transition at T s > T c . Fitting the Bg suscepti- 
bility data near T c to ([3]) with both y t and yh free, we find 
a second Bg transition at T = 0.4316(4), with exponents 
y L = 1.00(1) and y h = 1.94(1). The estimate y t = 1.00(1) 
agrees perfectly with the exact Ising value yt = 1, and 
the transition temperature is also entirely consistent with 
our estimate of the chiral transition temperature T c . We 
therefore conjecture that at the chiral transition temper- 
ature T c , the spins on the Bg and Cg sublattices undergo 
an Ising transition separating two QLR-ordered phases. 
The value of yh = 1.94(1) obtained from the fits of Xb 8 
near T c is clearly not equal to the usual Ising value of yh 
however; we therefore further conjecture that this value 
is a newly observed critical exponent of the Ising univer- 
sality class. Determining the exact value of yh = 1.94(1) 
remains to be further explored; we note however that 
from the Kac table, it may be ^ sa 1.9479. . .. For il- 
lustration, the inset of Fig. |4] shows xb s L~ 91 ^ 4s vs T. 

As further evidence that the Bg sublattice undergoes 
transitions at both T s and T c , Fig. [5] shows the specific- 
heat-like quantity C_b 8 . We clearly observe a peak near T s 
and a divergence near T c as L — > 00. For comparison, the 
inset shows Ca 4 , which does not show any peak near T s . 
To illustrate that the phases of the Bg spins on either side 
of T c are indeed both QLR ordered, we display in Fig. [6] a 
log-log plot of the xb b L~ 2 data vs L^ 1 for various values 
of T. We clearly observe an algebraic divergence of xb s , 
with a T-dependent exponent. 



Bisected-hexagonal lattice. — Analyzing Xc and Q c in 
the same manner as for the triangular and UJ lattices, 
we find that the chiral order parameter undergoes an 
Ising transition at T c = 0.39137(8). Similarly, an analy- 
sis of the sublattice susceptibilities shows that the spins 
on each sublattice undergo a transition from disorder to 
QLR order at a common value of T s = 0.747(2). The 
magnetic transitions are again consistent with the KT 
universality class. As for the Bg and Cg sublattices on 
the UJ lattice, xa 4 and xb b show that the A4 and Bq 
sublattices of the BH lattice also undergo an Ising transi- 
tion at T = T c . Near T c , the specific-heat-likc quantities 
Ca 4 and Cb 6 are observed to be divergent for L — > 00, 
while Cc 12 seems to be a smooth function. By contrast, 
near T s , a peak occurs in Cc 12 which is absent in Ca 4 
and Cb 6 ', see Fig. [7] These observations suggest that the 
magnetic transition on C12 is a generic KT transition, 
and that the QLR magnetic order on A4 and Bq in the 
region T c < T < T s is likely induced by the long-range 
correlation length on C\2- 

Discussion. — In this Letter, we have studied the XY 
antifcrromagnet on three plane triangulations: the tri- 
angular, UJ and BH lattices. Each of these lattices is 
tripartite and Eulerian (all vertex degrees are even). In 
each case, we find that the chiral order parameter un- 
dergoes a standard Ising order/disorder transition. This 
behavior should be generic on all Eulerian tripartite tri- 
angulations of the plane, since the set of ground states 
will share the same Z2 chiral degeneracy. However, as 
we have shown, the nature of the magnetic transitions 
appears to be strongly dependent on the specific lattice 
topology. 

In addition to the intrinsic theoretical interest of these 
results, we note the AFXY model on the UJ lattice could 
be realized experimentally using Josephson Junction Ar- 
rays in an entirely similar manner to the triangular- 
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FIG. 7: (Color online) C Ai versus T on the BH lattice. The 
inset is for Cc 12 - 



lattice model. Therefore, the results that we have out- 
lined above should be experimentally observable. 
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